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Noise Generation in Hot Jets 


Abbas Khavaran 
ASRC Aerospace Corporation 
Cleveland, Ohio 44135 


Donald C. Kenzakowski 

Combustion Research & Flow Technology, Inc. (CRAFT Tech) 
Pipersville, Pennsylvania 18947 


A prediction method based on the generalized acoustic analogy is presented, and used to 
evaluate aerodynamic noise radiated from high speed hot jets. The set of Euler equations 
are split into their respective non-radiating and residual components. Under certain 
conditions, the residual equations are rearranged to form a wave equation. This equation 
consists of a third-order wave operator, plus a number of non-linear terms that are 
identified with the equivalent sources of sound and their statistical characteristics are 
modeled. A specialized RANS solver provides the base flow as well as turbulence quantities 
and temperature fluctuations that determine the source strength. The main objective here 
is to evaluate the relative contribution from various source elements to the far-field spectra 
and to show the significance of temperature fluctuations as a source of aerodynamic noise in 
hot jets. 


I. Introduction 

T FIBRE is ever-increasing demand on the aircraft manufacturers to achieve meaningful levels of noise reduction 
in order to meet the stringent regulations that are being set for noise emission. To address the noise concerns 
near airports, NASA introduced a stretched ‘■‘pillar’ goal in 1997, with a stated goal of reducing the perceived noise 
levels of the future aircraft by 20 EPNdB by 2022 relative to 1994 technology. Jet noise is of particular interest 
because of its dominant contribution to the community noise at the high thrust conditions during takeoff. As part of 
an ongoing research effort at the NASA Glenn Research Center, a statistical-based jet noise prediction methodology 
(JeNo code) is being developed to provide the state-of-the-art in the source modeling and propagation. Since jet 
mixing-noise is broadband, with frequencies that cover significant levels of noise within three octave band, a 
successful predictive capability needs to incorporate sufficient physics to produce a reasonably accurate spectrum 
over a broad range of conditions. The goal is to achieve practical levels of accuracy in noise prediction from this 
more generalized modeling methodology to help designers with concept evaluation and down selection within a 
reasonable time frame. 

Predicting and suppressing noise from heated jets is of particular interest in aeroacoustics since this is the 
condition under which real engines operate. In an earlier study [Ref. 1], an acoustic analogy methodology was 
presented based on the premise that an identified fluctuating velocity/enthalpy source term was a major heat-related 
source contribution in hot jets. The present paper aims to further examine individual source terms in an 
inhomogeneous variable density Pridmore-Brown wave equation for their relative contribution to the far- field jet 
noise spectra. Lighthill [Ref. 2] defined the equivalent sources of aerodynamic noise as the double divergence of 

the stress tensor T.. =pv i v j +(p — C x p)8 jj —T ij , where p denotes the density, v = (v l ,v 2 , V 3 ) is the fluid velocity, 
p is the pressure, C m is the ambient speed of sound, and r. . denotes the viscous stress tensor. This definition has 

been a subject of wide-spread scrutiny and interpretations and has led to scaling laws that relate the far-field sound 
to flow parameters such as velocity and temperature in a variety of forms. The viscous stresses are usually 

considered to be unimportant in noise generation and the momentum flux term d 2 pvv Jd 2 x j x j is regarded as the 

primary source of sound in isothermal jets. Heat addition has a multitude of effects on noise that depends on the jet 
velocity, frequency, and angle. Experimental observations [Refs. 3,4,5] show that at high speeds (i.e., above 
acoustic Mach numbers of 0.9), heat addition results in noise reduction at all frequencies. At low speeds, nominally 
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acoustic Mach number of 0.50 and below, it amplifies the low to mid frequency jet noise with minimal effect at the 
high frequency. In between, e.g. at acoustic Mach number of 0.70, the spectrum shows a cross-over relative to the 
unheated case that indicates an enhancement of the low-frequency noise and slight weakening of the high-frequency 
amplitude. 

These observations have compelled many researchers in the field to examine the pressure/density 
difference (p —C^p) that appears in Lighthill’s stress tensor [Ref. 2], as a potential contributing source and 

formulate scaling laws that resemble the usual (7 8 law suggested based on the first source term. Some recent 
examination of jet noise data [6, 7] conclude a more general power law AU " , where the amplitude A and exponent 
n are parameters that depend on angle and jet stagnation temperature. Morfey et al. [8] revised Lush and Fisher’s 

[9] version of entropy source term and offered a new two-scaling hypothesis, in favor of U 6 and (7 4 power laws, 
and with some dependence on the mean temperature gradient. Additionally, they attempted to account for the mean 
flow effects by using a simplified form of the Geometric Acoustic approximation. Lilley [10] proposed a second 

order wave equation for the pressure fluctuations, similar to Lighthill’s, but replaced ( p —c\p) with a different 
expression that explicitly displays its isentropic and non-isentropic components. His new contribution to the total 

acoustic power consisted of a Dipole term with U ' dependency, which also multiplied the variance of the total 
enthalpy. He provided some estimate of the relevance of the new source term without actually calculating the 
required enthalpy fluctuations 

The acoustic equations described in the following section are derived from equations of motion after each 
flow variable is decomposed into its base flow and radiating components. The Euler equations, rewritten within the 
framework of Goldstein’s generalized acoustic analogy [Ref. 1 1], are expressed as two set of equations that govern a 
non-radiating background flow plus its residuals. The residuals components are simply the difference between the 
original dependent variables and their base-flow equivalents. These are later rearranged to form alternate analogies, 
and equivalent sources of sound are identified with the non-linear terms in each analogy. Of key significance in the 
current study, all analogies include, among other factors, a fluctuating velocity/enthalpy term as a heat-related 
source component. 

The non-radiating base flow alluded to above is taken to be the jet mean flow as calculated from a 
specialized RANS solver that predicts the variance in total temperature (or enthalpy) in addition to the standard 
turbulence-related parameters. The approach, which is discussed in [Ref. 12], makes use of a baseline k-e turbulence 
model demonstrated in past work to predict correct mean flow velocity mixing and turbulent kinetic energy for 
heated and unheated subsonic round jets. A generalized two-equation scalar variance model is available [Ref. 12] 
within this turbulence model framework, which utilizes an independent dissipation rate equation as well as the 
locally available turbulent velocity and time-scale information. Determination of the scalar variance variable is 
principally achieved by selecting the appropriate mean flow quantity gradient for its production source term 
(Appendix-A). In the past, this model has been used to predict the local variations in the turbulent Prandtl number 
and to study its impact on jet thermal mixing using an energy variance approach that tracked the variance in static 
temperature. The current study traces the variance of the total enthalpy instead as the parameter of interest for 
subsequent jet noise predictions. The stagnation temperature fluctuations are then deduced from this quantity. The 
difference between utilizing static versus stagnation fluctuation values for noise prediction becomes increasingly 
obvious with jet speed. Specifically, a noticeable static temperature variance appears within the jet shear layer for 
unheated conditions at high exhaust speeds, compared to a practically insignificant total temperature variance. From 
the standpoint of turbulence model generality, it is interesting to note that, at least for the jet cases studied, the use of 
total enthalpy variance does not significantly alter predictions for variable turbulent Prandtl number. The local 
turbulent thermal diffiisivity is partially determined from the thermal time-scale; the predicted values for the 
variance dissipation rate are also affected proportionately by the scalar variance production term selection. 

The remainder of this paper is organized as follows. The governing acoustic equations and the sources of 
jet noise in the specialized case of a locally parallel mean flow are presented in section 2. A summary of the 
baseline scalar variance model governing the transport of total temperature variance and its dissipation rate is 
provided in Appendix A. Details regarding the model development can be found in [Refs. 12 and 13]. Section 3 
illustrates the sensitivity of stagnation temperature variance to heated and exit velocity conditions for a number of 
round jets. A comprehensive study of the relative strength of various source components and their far-field noise at 
90° is presented in section 4. The comparisons cover a host of subsonic jets in an attempt to identify the most 
significant sound sources for a wide range of conditions. Source statistics and modeling features appear in 
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Appendices B and C. The current prediction methodology is calibrated around two source components only, and 
sample results are presented in section 4. 


II. Acoustic Analogies 

The details provided here are for completeness and for drawing attention to the underlying assumptions and 
noise sources in each analogy. Spectral predictions use a third-order inhomogeneous wave equation. This equation 
assumes a locally parallel mean flow, and consists of a linear wave operator that accounts for the refraction of sound 
due to flow in-homogeneities, plus a host of source terms that are all non-linear in fluctuating variables (i.e., velocity 
and total enthalpy fluctuations). 

We neglect the non-isentropic effects considering that fluid viscosity and heat conduction play a relatively 
insignificant role in noise generation. The conservative form of Euler equations 


dp d 

— + pv, =0, 

dt dx. 

(la) 

d d dp 

— pv i + pv.v + — = 0 , 

dt dx dx . 

J 1 

(lb) 

3 3 

-(ph,-p) + —pv j h t =0, 

(lc) 


are linearized about a suitable background flow, and expressed as two set of equations that govern a non-radiating 
base flow and its residuals. The residual equations are the difference between the original Euler equations and the 
base-flow relations, and govern the small fluctuations that relate to the generation of sound. 

Assuming that the equation of state is governed by the ideal gas law 


p = fMT , 

each flow variable is divided into its mean and fluctuating components as 
P = P + P ', P=P+P, v. =v. +v', h = h + h' . 


(2) 


( 3 ) 


The over-bar is a time-averaged value, and tilde denotes a Favre-averaged (i.e., mass-averaged) quantity for all 
variables including the stagnation enthalpy 


1 , _ — 
h = h + — v , h - h + h ' , h' - 0 . 
2 


(4a) 


The following expressions are concluded for the Favre-averaged stagnation enthalpy h and fluctuations h' 


7 1 , ~ 2 , ’ 2 7 ' 7 ' , ' 2 1 '2 , ~ 7 

n t =h + -v +-v , h l =h+-v — v +v ( .v ( . . 
2 2 2 2 


(4b) 


The averaged static enthalpy satisfies the averaged ideal gas law 

h = cT, p = SRp f , (4c) 

and upon subtracting (4c) from Eq. (2), the residual part of the ideal gas law becomes 


NASA/CR— 2007-2 14924 


3 



(4d) 


3? 

p = — (ph' + p'h), h' = cT', 

c 

p 

where c p and c . denote the specific gas constants (i.e., C = V C v ). Equation (4d) shows that density and 

temperature fluctuations act in concert to generate pressure fluctuations. This is also seen from the equation of state 
by taking its total differential 


dp dT dp 
p T p 


(4e) 


If the gas were to expand slowly enough to make it possible for the temperature to remain constant, then dT 
would have to remain as zero, and according to (4e) pressure and density would have to relate as dpi p = dpi p . 
However, the expansion and contraction of sound waves happen on rapid scale. It is more appropriate to assume 
that the thermal energy has very little time to flow, and so the expansion occurs adiabatically. Additionally, when 
the viscous effects are small, sound generation should be regarded as an isentropic process, resulting in 
pressure/density relationship dp Ip — y(d pip). Therefore, in view of (4e), changes in pressure and temperature 
relate as dplp = y/(y - 1 )(dT/T) . 


The analysis shows that density fluctuations may also be replaced with their equivalent in terms of velocity 
fluctuations from the continuity equation. With the aid of momentum and energy equations, the pressure fluctuations 
are subsequently expressed in terms of velocity and thermal fluctuations. 

The set of Euler’s equations are averaged using the above definitions. The averaged equations, herein referred to as 
the base-flow equations are 


D oP = 0 , 

(5a) 

, dp d 

d Ap v ,) + — = "— pv,Vj, 

OX. ox . 

1 J 

(5b) 

dp d 

D ( ph ) = pv h . 

‘ dt dx 1 ' 

(5c) 


To obtain a self-contained system of equations for the base flow, energy equation (5c) is rearranged as 


D 


p(h + — v ) 
2 


dp d 777 1 — 

= pv .h — D (pv .v .) 

dt &. 7 ' 2 ° ^ 1 j} 


(5d) 


The base flow is solved by modeling the stress terms on the RHS of (5b) and (5d) and using the ideal gas law 
(y — Y)ph = yp. In a RANS-type mean flow we require dp I dt — 0 . Upon subtracting the above equations from 
the set of equations (1), the following relations (referred to as the residual equations) are found for the fluctuating 
variables 


dm . 

D p '+ — ^ = 0, 

0 dx. 

dp ' d _ de .. 

D m. + D (p'v.) + + ( m v, ) = — 

O l O l ' /-s /-n V J l J ’ 

OX. ox. ox. 

1 J J 

,8 , dp' d — y—. 

D o (P K + pK) + — (pvjh,)- — = —pv i h l . 

OX. ot ox 

J J 


(6a) 

(6b) 

(6c) 
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We have used the following definitions 


df d 

(7a) 

m . = pv. 

(7b) 

e ij = -{pv'y’j - pv’y’j ) . 

(7c) 


The energy equation is now written in a different form that is more closely related to pressure fluctuations p . Using 
the ideal gas law (4d), the first term on the LHS of the energy equation (6c) is recognized as 


( p'h, + pti) = P + mv + - (p v 2 - e,, ) . 

y - 1 2 


(7d) 


For later derivation of the convective wave equation, it is also more appropriate to express time derivative dp’ /dt 
in (6c) in terms of the operator D 


dp' , , dv _ dp' 

D p - p v. 

dt dx. dx. 


( 8 ) 


We substitute for dp' 7 fix,- from (6b) into (8), and subsequently place dp' / dt from (8) into (6c) to obtain the 
following energy equation 

1 d dv dv. i de,, Q 

D p -i (m h) + p — - + v m — - = —De jj —v. — (pv'.h', — pv'.h', ) (9a) 

y - 1 ° dxj J d Xj ' y dx, 2 ° 11 1 dx, dx . 1 ’ 

where 

h' t = h' + —v' 2 (9b) 

2 


The energy equation (9a) may also be written in a slightly different form using a new variable h' that is related to 
the actual fluctuation in stagnation enthalpy h', as 


t f j r , 1 r2 i i 1 r2 ~ i 

h = h =h + — v — v.v. 


Substituting (9c) into (9a), the energy equations becomes 


(9c) 


1 d dv. dv. i dv. d 

D p + ( m h) + p — L + vm — - --D e,, + e„ — (pv'.h' - pv'.h') (9d) 

y-1 ° P dx. y 1 ’ ^ dx. ' J dx. 2 ° 11 u dx. dx ^ 1 ° ^ J o) 

Alternate definitions for the dependent variables are possible at the cost of re-defining the source. For example, the 
first term on the RFIS of the energy equation (9d) may be combined with the leading term on the LF1S to define a 
new pressure variable [Ref. 1 1] as // = p ~(j- d)e kk H 
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dv. 


dv. 


7 - 1 


dv. d 


D p ( m h) + p — - + v.m . — - = (e.. e tA ) — (pv' h' — pv'.ti) . 

r - 1 ° fie. ' ' e dx. ,J 8x. y 2 kk,jJ dx. fit. 7 ° p 1 o) 

• i li il 


With respect to the generalized pressure p ' , the momentum equation (6b) becomes 


D m. + D {p'v.) + + -^—{m v ) = -^—{e —— e kk S.). 

0 ' ° lJ dx. dx J 1 dx 11 ' kk v 


y-\ 

2 


A. Special Cases of the Acoustic Equation 

Consider the case when the mean flow is parallel and a function of the span-wise coordinates only 
v - 8 { U (x 2 ,x 3 ), p - p(x 2 ,x 3 ), p - constant , 
subsequently 

d 8 

D =D = vU — . 

dt dx t 

Using the ideal gas law, we have 
{y - \)p h = yp = pc 2 . 

And the base flow momentum equation (5b) shows that 


pv .v = 0 . 

dx. ' ' 


Then according to Eq. (4b) 


h = h(x 2 ,x 3 ) + — U 2 (x 2 ,x } ) + - v' 2 . 

2 2 


And the base flow energy equation shows that 


— U — {pv 2 ) + -^— pv'.h' = 0 . 

2 fix, dXj 1 

Using the above results into the residual momentum and energy equations we conclude 

dp' „ dU de 

+ Dm + S .m . = — - , 

OX. ox. ox 

1 J J 

d ~ 

Dp '+ — {c~m ) = q, 

dx. 


where 


f 


- D (P V ’ 2 ) + (pv'jK )~ u y— ( pv[v'j ) 

2 dx. dx. 

J J 




q = -(r- 1 ) 

Source expression q may also be expressed in terms of h' as defined in (9c) 


(9e) 

(9f) 

(10) 

(11a) 

(lib) 

(12a) 

(12b) 

(12c) 

(13) 

(14) 

(15a) 
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q = -(,/- 1) 


( 1 n . /2 \ d , , dU^ 

— D(pv ) +—{pv j h o ) + {pv x v j )— 

y 2 OXj V X iJ 


(15b) 


Residual parts of the continuity (6a), momentum (13), and energy equation (14) are now rearranged to construct 
wave equations for either pressure or density. 

A form of Phillips’ equation is derived after eliminating factor m . between momentum and energy equations by 
taking the convective derivative of the latter and subtracting that from the gradient of the former 

8p ')_ 8 (7 1 8)6,1 

dx . 8x. dx . 

J 1 J 

Lilley’s second-order wave equation for acoustic pressure in a stationary’ mean flow [Refs. 10, 14] is deduced from 
Eq. (16) if we let U — 0 , c~ = c 2 m , and p — (note that fi\ — h[ in a stationary mean flow) 


) + 2c‘ 


8U dm, 
dx. dx, 


- + Dq. 


(16) 


D 2 p' (c 2 

dx. 


d 2 p' d 2 p' _d 2 (pvV.) y-\d 2 {p V ' 2 ) d 2 pv ' h’° 

c 2 Jdt 2 dx j dx j dx.dx j 2 c 2 dt 2 dxpt h m 


(17) 


where (y — 1 )h t = c\ . Since the sound speed was set equal to that measured at the far-field observer point, the 

mean density needs to be equal to the ambient density p — p_ = y pjc " as well. Therefore the actual jet density 
does not have a place in the source. 

Similarly, Lighthill’s equation for density fluctuations in a quiescent medium ( U = 0, c~ = c‘ , 
and p — p r ) is written after eliminating d~m Jdtdx . between continuity and momentum Eqs (6a) and (13) 


^2 f ^2 f 

dp 2 d p 

dt 2 “ dx dx . 

J J 



(18) 


Again p should be set equal to the ambient density p rl . Notice that unlike the conventional acoustic analogy, the 

base-flow equation is already subtracted from Eq. (18), and the sources on the right hand side consist of fluctuating 
terms only. The pressure/density difference that appears as a source in Lighthill’s equation is second-order in 
fluctuating variables, and may be written after eliminating dm ./dx x between continuity and energy equations (6a) 

and (14). A form of Lighthill’s acoustic analogy that accounts for the energy equation would be a third-order wave 
equation 


d_ 

dt 


d 2 P 2 

PP~ C ’ 


ay 

dx .dx . 

J J 


= P, 


dx.dx . 

l J 


— v.v. -(y- 1) 

dt 1 J ’ 


d r2 d r ) 
-i v,h i 


2 dt 


dx,. 


J 


(19) 
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A comparison of equations (17) and (19) clearly shows that pressure fluctuations may be viewed as a preferred 
parameter in writing the wave equation. 

The source components in equation (17) when calculated on a proper surface (i.e., U— 0) may be used in 
conjunction with the Kirchhoff approach or Ffowcs- Wiliams/Haw kings (FWH) surface integral methods to calculate 
the far- field acoustic pressure. The difference between the above source and the usual acoustic analogy (18) is that 
Eq. (17) requires time history of the velocity and the stagnation enthalpy on the surface. 

Rather than assuming a quiescent acoustic medium, it is more appropriate to include the mean flow effect by moving 
the linear component from the second source term on the right hand side of Phillips’ equation (16) to its operator 
side. To realize this, form the axial gradient 8/dx of the momentum equation (13), multiply that by dU/dx , and 
subtract the result from the convective derivative of Eq. (16) to obtain the third-order convective wave equation 


2 , 5 ~ dp' 

D D-p (c — ) 

dx. dx. 

J J 


~\~ 2c 


dU d 2 p 
dx. dx.dx. 

I 1 / 


de.. 


-D — (c — — ) + 2c 
dx. dx . 


dU d-e v 

dx. dx.dx. 
j i ' 


+ D 2 q. 


(20) 


Stress terms £?.. and q were defined in Eqs. (7c) and (15). This result, also referred to as an inhomogeneous 
Pridmore-Brown equation, is expressed in a more transparent form as 

Lp= r, (21) 


where source expression T was defined as the REIS of Eq. (20) and the operator L is 


, d ~ d ~2 

L = D (D 2 c 2 ) + 2c 2 

dXj dXj 


dU d 2 

dXj dx x dXj 


(22) 


Equation (20) may also be derived directly by rearranging equation (3.5) of [Ref. 11]. 

A comparison of equations (20) and (17) shows that the presence of a mean flow not only refracts the sound after it 
was generated, but it also modifies the sources of sound as evidenced by the REIS of these two equations. 


B. Discussion 

The acoustic equations described above are all consistent in their description of the aerodynamic noise sources. 
With the exception of Phillips’ equation (16) which still holds a linear term on the right hand side, the subsequent 
wave equations (17), (19) and (20) are summarized as Ltj=i $ , where the linear operator L acts on either pressure or 

density fluctuations, and the source s’ consists of components that are, at least, second-order in fluctuating variables. 
These equations were derived from the residual continuity, momentum and energy equations (6a, 13, 14); and as 
such their respective sources should be viewed as the difference between the more traditional analogies and their 
non-radiating base flow. In conventional analogies, source consists of both linear and non-linear components 
[Refs. 15,16,17]. In the following, additional simplification is provided prior to an assessment of individual source 
elements. 

C. Further Simplification of the Acoustic Equation 

The source expression in the convective wave equation (20) is further simplified once we insert the 
assumption of a locally parallel mean flow into the base-flow equations and conclude that Reynolds stress pv'.v' has 
a Favre -averaged value that is independent of the source location. While this simplifications leaves out the gradients 
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of pv'.v from the source elements in Eq. (20), it is still realized that in a real spreading jet, base-flow momentum 

and energy equations are both inhomogeneous, which implies that the true sources are the difference between the 
fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. Subsequently we let 


de.. 
']_ 

dx. 

J 



(23) 


In a RANS-based jet noise prediction, the base flow is usually the averaged mean flow. The two-point statistics 
between various fluctuating variables are modeled and calculated from their respective variance combined with the 
time- and length-scales of the local turbulence. Thus the density factor p appearing in the numerator of each 

source component in Eq. (20) is replaced with its mean value p , which practically amounts to neglecting terms that 
are higher than second-order in fluctuating variables. Dividing equation (20) by constant y p one finds 


Ln = T, 


n 


YP 


where the equivalent sources of aerodynamic sound are 


( j_ 

2 r t 


Y = D 


ii 





dx.dx . 

* J 


V 


f I 

_~8U_ 

^2 r f 
d V iVj 

dx. 

J 

dx t dx. 


>~\2 t ! . 

0v ' v '-i 

2 h 


self 


J 


ii 


hi 


t t 

, VjV, 
I 2 / 1 J 


1 dp d 


V- ( x )+< ^V <v:v;) 


-D 


ii 




dx : 


h dx. p dx i 


shear 


enthalpy 


(24) 


(25) 


Three pairs of brackets in equation (25) designate three distinct sources. The former terms in each bracket are 
identified as self, shear, and enthalpy sources respectively. Self and shear source terms are second-order in 
fluctuating velocity, and are each complemented by a second term involving the mean flow enthalpy which is not 
zero even when the flow is isothermal, but is shown to be small relative to the former. Source terms that consist of 
the mean density gradient are negligible in unheated jets, but could potentially become important as jets get hot. 

The second source component within the enthalpy noise group was originally part of the first source term in 
equation (20); and is now grouped with the heat-related sources. A detailed comparison of the individual source 
terms of Eq. (25) is presented in the next section. It turns out that the major difference between the unheated and 
heated jets is due to the fluctuations in stagnation enthalpy that appears in the third bracket above. This term is first 

order in velocity fluctuations, and its spectral shape behaves as in the low frequency limit, and its high- 
frequency roll-off could be modeled to behave as a> , similar to the self noise. 


III. Subsonic Flowfield Simulations 

A total of eight subsonic cases within the Tanna matrix [Ref. 18] were considered for simulation of a 2-inch 
diameter jet using the total enthalpy variance model (Appendix A). The nozzle exhaust conditions as defined by set 
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points within the Tanna matrix are summarized in Table I. In past work, these cases were studied extensively [Ref. 
19] to better understand and predict turbulence quantities in the developing jet region, where most of the jet noise is 
produced. Comparisons were made with available PIV data collected at the SHJAR facility [Refs. 4,20] at 
NASA/Glenn Research Center. Sample results that compares energy variance model with the total enthalpy variance 
model were presented in a previous study [Ref. 1]. 

Simulations use an axisymmetric 316x171 computational grid, which includes the internal nozzle region 
well upstream of the exit plane. The computational domain extends 100 jet exit radii in the streamwise direction and 
25 jet exit radii radially, sufficiently distant to minimize boundary condition placement impact on the shear layer 
entrainment path-lines. Nozzle flow boundary conditions were prescribed uniformly as inflow stagnation conditions 
for pressure and temperature imposed well upstream of the nozzle exit. Solution grid resolution sensitivity and 
convergence criteria were discussed in [Ref. 1]. The same reference also compares fluctuations in total and static 
temperatures and shows that the former parameter remains consistently small in unheated jets while the latter could 
grow to become significant with jet speed. The predicted centerline decay of mean axial velocity and turbulent 
kinetic energy are plotted in Figure 1 as a function of Witze axial parameter x for subsonic jets [Ref. 21]. The 

length of the potential core collapses at xj r. = 1 in agreement with the theory - and the peak TKE on the centerline 

is at xj r. =1.5, and shows slight increase with jet temperature. 




w 



U. 

-0.16 — 

c. 




(%-°' 22 x 

Pi 


(26) 


Contour plots in Figures 2 and 3 show the variance in total temperature for the first 15-diameters of the jet as a 
function of jet speed and temperature. The peak temperature fluctuations increase with both speed and temperature. 
Similar calculations in unheated jets show that peak levels of 0.20x10-2 and 0.70x10-2 for cases sp03 and sp07 are 
relatively insignificant within the context of the heated jets. 


Table I. Tanna Matrix Set Point Conditions 


S.P. 

NPR 

T It 

t,J / °° 

t It 

J / °o 

M * 

ac 

M 

J 

03 

1.197 

1.000 

0.950 

0.50 

0.51 

23 

1.103 

1.810 

1.760 

0.50 

0.37 

42 

1.066 

2.750 

2.700 

0.50 

0.30 

07 

1.861 

1.000 

0.840 

0.90 

0.98 

27 

1.361 

1.922 

1.760 

0.90 

0.68 

46 

1.225 

2.861 

2.700 

0.90 

0.54 

29 

1.900 

2.114 

1.760 

1.33 

1.00 

49 

1.692 

3.138 

2.700 

1.48 

0.90 


*M ac =C7,/ C[ 
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Figure 1. Centerline velocity and turbulent kinetic energy vs. Witze axial parameter in subsonic jets. 



ARN - sp29 


Figure 2. Variance in total temperature normalized 
with exit total temperature at static temperature ratio 
of 1.76 and jet acoustic Mach number of: 

0.50 (top); 0.90 (middle); and 1.33 (bottom). 



ARN - sp49 


Figure 3. Variance in total temperature normalized 
with exit total temperature at static temperature ratio 
of 2.70 and jet acoustic Mach number of: 

0.50 (top); 0.90 (middle); and 1.48 (bottom). 
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IV. Acoustic Results 

In the first step, the relative strength of various source terms in the governing wave equation (Eq. 24) will 
be compared at jet exit conditions shown in Table I. This exercise helps to identify sources that are relatively more 
important in the noise generation process. The remaining sources will be considered as small and neglected in the 
final calibration of the prediction methodology. 

A. Source Strength Evaluation 

Jet noise spectra from various source components are compared at 90° (see Eq. C5, Appendix C). The source 
correlation coefficient / is replaced with 7 nu in an isotropic turbulence as described in Appendix B. Both non- 
compact (Eq. B9) and compact (Eq. B21) variations of the source are considered to in order to ensure that important 
noise components are not inadvertently overlooked as a result of the specifics of the modeling. 

Each row within the column matrix Eq. C5 represents, in an orderly fashion, a source element in Eq. 25. The first 

and third source terms within the shear noise group are identically zero at 9 = 90“ (i.e., rows 3 and 5 in Eq. C5); and 
the remaining five sources are integrated over the jet volume starting from its exit plane. The turbulence length- and 

time-scales are calculated as £ = Ck 5 / s and T = CkI s for all source elements. Figure 4 shows the relative 

significance of each source component when a non-compact source model is used. Comparisons point toward self2 
and shear 2 as sources that are relatively small under all conditions, whereas enthalpy 1 and enthalpy 2 intensify with 
heat and become serious contenders to the sefl under heated conditions. Between the two heat-related sources, the 

former scales with (o' at low frequency and dominates most of the spectrum - while the latter scales as Of and 
should be viewed a potential source only at the very early stage of the spectrum. Component predictions shown in 
Figure 5 use a compact source model (Eq. B21). A similar relative strength among sources is observed, although the 
spectra appear slightly less broad relative to their respective counterpart in Figure 4. 


B. Spectral Predictions 

Numerical results are shown for eight subsonic set points described in Table I, and are compared with the 
lossless, narrow-band measurements of Ref. [4] at a distance of R/D =100. In addition to the standard flow 

variables, which include the turbulent kinetic energy and time- and length-scales, this modified version of the flow 
solver permitted input of the local total temperature variance as predicted by the modified scalar variance model 
(Appendix A). Computations use only two source components denoted as selfl and enthalpyl, as represented by the 
first and the sixth rows in Eq. C5 (Appendix C) 


p 2 (x,y,a> ) 


p 2 J nu k 4 {l -M s cos 6f {cos 2 6 + QsmOf 
(AxR) 2 (1 -Af cos 0) 2 


A(cos 2 0 + QsiwO) 2 


+5(1 -M' cost?) 2 


16 k h 2 _ 


Xu +<u/'"r* 

m = 0 


(27) 

Outside the usual calibration constants C, and C r that determine the length- and time-scales, the final calibration is 

achieved by selecting a pair of constants A and B in Eq. (27), which practically amounts to using a linear 
combination of the two sources. This is necessary due to the difference in scales between velocity and thermal 
fluctuations as well as possible cross-correlation between various sources. Parameters A and B were determined by 
calibrating the 90° predictions across two set points SP07 and SP27 (Table I). 

Spectral results are shown in Figures 6, 7 and 8 for the unheated jets as well as jets at temperature ratio of 1.76 and 
2.70. Predictions exhibit a reasonably good agreement with data at angles not very close to the downstream jet axis. 
The near-axis results weaken gradually with increasing jet speed and/or its temperature. The deteriorating quality of 
shallow-angles predictions is more likely attributed to the neglect of jet spread in the Green’s function calculations 
which also leads to a dismissal of the causal Green’s function. Fine-tunings such as cross-correlation between 
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sources, or more elaborate statistical models with multiple length-scales and/or non-separable correlations are 
unlikely to correct the near-axis deficits observed in a parallel flow approximation. The causal solution to the 
homogeneous parallel-flow equation (20) is known to become unbounded downstream, and is usually discarded as a 
viable Green’s function. A second bounded Green’s function is found only when the mean flow supports a small jet 
spread [Ref. 22]. In that case, a complete solution to the sound field is written as a convolution of the source 
elements discussed here with both components of the Green’s function. 


V. Summary 

Identification of the true aerodynamic noise sources is an important aspect of the aeroacoustics research in 
order to help engine designers with their noise mitigation efforts. Prediction of jet noise from empirical models is a 
quick and practical way of noise assessment, but is also limited to the bounds of the data-base. The acoustic analogy 
derivation discussed in section II provides a coherent noise prediction framework that recognizes the sources as the 
difference between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. As shown through 
the component comparisons of section-IV, source elements identified as self ! and enthalpyl in Eq. (25) emerge as 
the major contributing components to the sound field. This study also highlights the benefits of providing the total 
temperature fluctuations as a formal input parameter for jet noise prediction from a RANS tlowfield simulation. 
The current approach in extending a baseline scalar variance model to furnish this additional input in a generalized 
manner is exploratory but appears promising. Needless to say, additional experimental measurements are required 
to verify the RANS model’s quantitative prediction for the temperature fluctuations across a range of propulsive hot 
jet conditions. Statistical models proposed in the source modeling, although intended to capture the main features of 
various sources for their strength and spectral roll-off, also require validation. It is hoped that further refinement of 
the modeling, in conjunction with the instability-related noise component, i.e., convolution of the causal Green’s 
function with the sources described here, could improve spectral predictions at angles close to the jet axis. 
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Figure 4. Effect of jet temperature on noise sources in Eq. 26. 

Predictions are at 90° observer angle and use a non-compact source model (Eq. B9). 
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Figure 5. Effect of jet temperature on noise sources in Eq. 26. 

Predictions are at 90° observer angle and use a compact source model (Eq. B21). 
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1 0Log (Uj D'j p 2 / ), dB 1 0Log (Uj D'J p 2 / P 2 ), dB lOLog (U, D ] p 2 / P 2 ), dB 






1 0Log (Uj D'j p 2 / P 2 ), dB 1 0Log (U : D'] p 2 / P 2 ), dB 





Appendix A - Enthalpy Variance and Dissipation Rate Model 

The high-Reynolds-number form of the equations governing the transport of the total enthalpy variance and its 
dissipation rate is based in concept on the work ofNagano and Kim [13] for generalized turbulent thermal 
diffusivity 


D 

~Dt 


pk t 



a + - 
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> kt J 



+ p k, - 2 PS, 


(Al) 
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ps, 
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ps 


— C 


d 4 


\p£ t 


(A2) 


Here, k and S denote the turbulent kinetic energy (TKE) and dissipation rate obtained from a companion two- 
equation turbulence model, the operator D /Dt is defined in a similar fashion as D in Eq. (7a), a and a t are the 

molecular and eddy diffusivities for heat, and Pi is the production term used in the TKE equation. Values utilized 
for the modeling constants in the variance dissipation rate equation are based on results described in [Refs. 1,12] and 
are listed in Table . 


Table Al. Enthalpy Variance Dissipation Rate Coefficient Constants 


c dl 

C d2 

c d3 

C M 

Gkt 

&£t 

2.0 

0.72 

2.2 

0.8 

1.5 

1.5 


Physical interpretation of the quantity k t tracked within the CFD simulation is based on selection of the flow variable 
gradient used for the turbulent variance production term P k For example, Nagano and Kim originally solved 

directly for the static temperature variance. Recent numerical studies modeling helium jets and afterburning rockets 
plumes indicate improved model generalization using an internal energy variance approach. Further extension of 
the framework to model the total enthalpy variance was found necessary to improve jet noise prediction. 

Two modeling approaches were studied for use in hot jet noise prediction [Ref. 1]. The first utilized the internal 
energy variance as described in [Ref. 12]. The second modeled the total enthalpy variance, with appropriate variable 
substitution for the variance dissipation rate equation. The production term utilized in the total enthalpy variance 
model is 


n = 2«, 


dh, dh, 

dXj dXj 


(A3) 


where e is the specific internal energy and h, is the specific total enthalpy. For a calorically and thermally perfect 
gas, the specific heat constants are invariant with temperature, and so the temperature variance can be expressed as f 


k , = KK = c p 2 (7777) , 


(A4) 


and c is the specific gas constant. When extending the variance model to predict thermal eddy diffusivity, an 

expression for the thermal turbulent time scale is needed. Based on results from past work, the following expression 
was used: 



V 


(C ^ 

fjfcV 

— max 


? 

— 

U J 



K C * J 

U J 


C M = 0.09 


(A5) 


The thermal eddy diffusivity, applied to the turbulent diffusion term in the mean energy equation and the variance 
production term, is written for model closure as 
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«/ = C xP kT , ■ 


(A6) 


As recommended by Chidambaram et al. [23], value of C A — 0. 14 was used. The local turbulent Prandtl number is 
then defined as the ratio of the eddy viscosity to the thermal eddy diffusivity 

Pr, =//,/«, . (A7) 


Appendix B - Source Model 

The far-field sound spiral density per unit volume of turbulence at y is 


p 2 (x,y,co) - |G(x,y,a>)|~ J J q(y, , r)e ,alT e ' k4 dzd% (Bl) 

I -o 

where q denotes a two-point space-time correlation between fluctuating variables at source points A and B separated 
by space d, and time T , and G is the relevant Green’s function [Ref. 24]. For simplicity, we consider the 
turbulence as isotropic as described by Batchelor [Ref. 25], and describe a two-point space-time correlation of the 
velocity components (notation U t is now used in place of v\ as a fluctuating velocity component.) 

(B2) 

fXM-jz, f 1 = £+£+£. 

Consider a fourth-order time-delayed space-time correlation per unit volume of turbulence at y 


R m (y^’ T ) = U i U j U k u 'e • ( B3 ) 

The spectral density associated with the axial component of the above tensor is the four-dimensional wave -number 
frequency spectrum function 

Inniy^l = JJXnO^I. r ) e“°‘ T e * 4 drd% . (B4) 

I r 

The above integral was written in a frame moving with source convection velocity U , and of is the source frequency 
in that frame and is related to the observer frequency through the Doppler factor. If turbulence correlation coefficients are 
quasi normal, then the axial component of the above tensor becomes u { u { u[u[ = 2 (u { u[) 2 and subsequently 


fun (y, (o s ) = 2(«r ) 2 J j [(/ + -z Zf ') - \ f f e ' a Te ' ki drd l ■ 


4 r 


2 # 


(B5) 


Similarly, consider the two point space time correlation for the enthalpy/velocity source term 


K 


h' 


1 


= ZjriU'hXW. , 

h h h 


(B6) 
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(B7) 


3 j(y,co s ) = J J 3..(v,|, T)e im ‘ r e ik4 dzd% . 

J r 

As suggested by measurements [Ref. 26], we assume that exponential functions govern spatial as well as temporal 
decay of velocity correlations 

/(£,r) = exp(-^-|r/r o |), A(r) = exp(-|z7r o |). (B8) 

Following Proudman [Ref 27] the correlation components are considered in the direction of the observer, i.e., the 
axial component of separation vector aligns with the direction of the wave number k (i.e., k .q — kq { ) . 

4f 3 ~ , r+a> r r 

H(g?) = J exp(— 2 I — I +ico s T)dT = - — / , (B9) 

5 n r o l + (®r o /2) 


where the non-compactness factor N t will be defined subsequently. Using Millionshtchikov’s hypothesis [Ref. 25] 
and noting that the autocorrelation of (uh' ) is zero in isotropic turbulence. 


( U ih'o)A(Uj h 'o)B =(Ui) A (Uj) B (K)AK) b ■ 

Now if the space/time correlation governing the enthalpy fluctuations is represented as 

then from Eqs. B6, B7, and BIO we find 
Sj j (y, (o s ) = ^-u* {{[(/ + 

n | r 


’) - YJ- / ’] g ($, r)e^e^dTdl , 


(BIO) 


(Bll) 


(B12) 


Experimental measurements show that both heated and unheated jets exhibit similar high-frequency decay. Ideally 
it is desirable for the correlation coefficients B4 and B7 to provide similar asymptotic behavior. In the absence of 
any experimental data, as a first approximation we assume that enthalpy correlation function g(q, r) decays in a 
similar fashion as f(iq, r) . From Eq. (B12) it is shown that 


" 3 - «Y 


(y, a ) = — - (m, ) — H (a s )N 2 (kl), g(£ r) = /(£ r ) . 


2 n* ' h 

The non-compactness factors 

5 


N l (2nx) = 


8/ 5 


3 tan - (x) - x 


5X-+3 

0+z 2 ) 2 


’ N 2( 2 *X)= 1 2 2 , 

(1 + / ) 


decay as X 5 an d X 4 respectively as / — >■ go . 

The spectral shape function due to either self or enthalpy source components is obtained from 


(B13) 


(B14) 
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(B15) 


F(co) = co 4 H(co)N(— t ) , 

which shows a low-frequency roll-off F(a> ) oc cq 4 as co —> 0 . The high-frequency decay is either F(o)) oc oj 3 or 
CO 2 depending upon factor /V [ or N 2 (see Eq. B14). The spectral density from above two sources relate as 


5 N 2 (kQ ]_h' 0 J 
8 N t (k£) u f /? 1111 


The ratio N^(kl)/ N^kl) remains as 1 .0 for most of the spectral range and somewhat beyond the peak, with the 

exception of the tail-end of the spectrum (Figure Bl). In evaluating the relative source strength, this difference has 
been neglected and the correct high-frequency behavior has been selected as 


5 

8 uf h 2 ' 


(B17) 



Figure Bl. Non-compactness factor for self and enthalpy source terms (Eq. B14); 
solid line, N t ; dashed-line /V, . 

In the following discussion, it is shown that alternative function forms are available that are compatible with the 
exponential functions, and provide similar high-frequency decay for both sources discussed above. 


Alternative Source Modeling 

A less complicated way of describing the statistical properties of turbulence assumes a compact eddy approximation 
in the limit of zero wave-number (i.e., k — > 0 ), in which case the non-compactness factors defined in Eq. B 14 
become identically equal to 1.0. If the correlations were expressed as separable, the spectral density would strongly 
depend on the temporal function of the correlation for its high-frequency roll-off, and very weakly on its spatial 
form, as the spatial Fourier transform converts into volume integration. Subsequently, similar spectral 
characteristics would be expected for both sources. 

For example, one may consider a temporal function such as (see Figure B2) 


h(t,V,a ) 


(at) v K v (at) 

2 v - l T(v) ’ 


h{t , \,a) = e at 


(BIB) 
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where time t is replaced with its absolute value when it is negative, and K y denotes a modified Bessel function. 
This function is void of a cusp at its origin and its Fourier transform scales as 0 ) 1 at high frequency 


rr , . - r°° . . . , ^V^r(v' + |) 1 

h (co) = 2\ h(t,v,a)cos{cot)dt = 2 — 7 -- 2/ 2 , r+/ • (B19) 

Jo a T(v) (1 + a 2 /a 2 y A 


The spectral shape function F(co) — O) H(co) which would be similar for both correlations B5 and B12, is now 
modeled by selecting parameter v . With v — 3 it is shown that 
F(co)^co a , asty— >0 

3 (B20) 

F(a>) — > a> , as co — > oo 



Figure B2. Function h(t, V, a) as defined in Eq. (B18); dashed line ( V — 3, a — 6 ); solid line exp(— 2t) . 


When v - 1 1 / 4 , the high-frequency roll-off follows F(co) — > CO /2 , and source correlation coefficients of interest 
become 


aF — 

Am =— r(«i 2 ) 2 #(®’)> 
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Appendix C - Source Green’s Function Convolution 

A stationary point source with frequency (O and location X s (superscript s denotes a source location) is considered 
in defining the Green’s Function (GF) 

/ - r 2 p- il0, fi( x - x s 1 


and 

1 . x 

G (x,x s ,co) = e ' k(R ~ x ' ms0) ^ f m (r , k, 9 ) cos m(cp - (p') , (Cl) 

4 ncoR m= Q 

where f m (r s , k, 9 ) is solved numerically as a solution to second-order compressible Rayleigh operator [Ref. 24], 
corresponding to mode number m, wave number k = col (to avoid confusion, notation K will be used to denotes 
TKE in this appendix), and observer angle 6 
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dr 

o 2 dr 
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, , , m ~ 

k~(Q>~ — cos 0 ) — 


m r\ 

g = 0 , 


O 2 = — (1 -M s cos 0 ) 2 , g m = (1 -M 5 cos O ) 2 f" . 

Poo 

The relevant GF associated with a moving singularity with source frequency (0 s and convection velocity iU is 


L{Ge iM ) = — 
Dt L 


ci^Six, -U c t)S(x t -x t ) 


, -* (1 - M s cos 9) rm , , , ^ 

G(x,x ,co) = e >/ (r ,k, 0 )cosm((p-(p ), 

4 nR (1 - M c cos 0 ) 


(C 2 ) 


where M s =U(r s )/c and M —U /c are the local acoustic Mach number and convection Mach number 

respectively. We make the simplifying assumption that the source components on interest are uncorrelated. As 
usual, the spatial derivatives are moved from source to the Green’s function. For example 

f ii ii . f (i 

(i — > U U . . 

dXjdXj ' 1 dx i dx j 


Some level of simplification is achieved when we consider components of turbulence fluctuations in the direction of 
the observer 6 as suggested by Proudman [Ref. 27]. As an example the first source elements within the self and 


enthalpy category (see Eq. 27) are designated as D 


8 ~u u d h' 

and D~ ( u —G) where 77 is the component of the 

dT]d?i dr] n h 


spatial separation of the correlation in the direction of the observer (see Figure). 
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The required spatial derivatives of the Green’s function are 

Q = cos 7 Q . + sin 7 Q 

,77 4 

Q = cos' 76,, + sin 2 76 + 2 sin 7 cos 76. 

>77 ,11 ,rr ,1 r 

where 6 , = —ik cos 7 6 and the radial derivative 6 r applies to the function f" (r, k, 7 ) . 

The following simplifications are introduced into the radial derivatives using the high-frequency solution to the GF 


fr = ~ikQf " , f:=-(kQ) 2 f m , Qr (1-M'cos7) 2 -cos 2 7, 

P . 0 


(C3) 


A more accurate representation of function Q should account for mode number 


Ql = — (1 -M s cost?) 2 - cos 2 7 -( — ) 2 . 

kr 


P '« 


(C4) 


Since the zero' 1 mode usually makes the major contribution to the GF, the last term on the RFIS was neglected in 
approximating the derivatives of the GF. Subsequently, the following expression are derived for the far-field 
spectral density per unit volume ring-source at y 


P 2 (x, y, co) = 
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(C5) 

Each of the seven components shown in the above column matrix represents the contribution to the sound field from 
a source of Eq. (27) in a similar order. The source correlation coefficient / ^ is the same as 7 im in either non- 

compact (B9) or compact-eddy (B21) representation, and K — 1.5m 2 is the turbulent kinetic energy. The turbulence 
length- and time-scales are calculated the usual way from the turbulent kinetic energy and its dissipation rate as 
t - C, k 1 ' 5 / £ and T = Ck/s. 
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